Persistent current noise and electron-electron interactions 
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We analyze fluctuations of persistent current (PC) produced by a charged quantum particle 
moving in a ring and interacting with a dissipative environment formed by diffusive electron gas. 
We demonstrate that in the presence of interactions such PC fluctuations persist down to zero 
temperature. In the case of weak interactions and/or sufficiently small values of the ring radius 
R PC noise remains coherent and can be tuned by external magnetic flux piercing the ring. 
In the opposite limit of strong interactions and/or large values of R fluctuations in the electronic 
bath strongly suppress quantum coherence of the particle down to T = and induce incoherent 
$a;-independent current noise in the ring which persists even at = when the average PC is 
absent. 
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I. INTRODUCTION 

Persistent currents (PC) in normal meso- and nanor- 
ings pierced by external magnetic flux represent one of 
the fundamental consequences of quantum coherence of 
electron wave functions which at sufficiently low temper- 
atures may persist at distances comparable with the sys- 
tem size. While the average value of PC was intensively 
investigated both theoreticallj*i and experimental!}^ dur- 
ing last decades, little was known about equilibrium fluc- 
tuations of PC. 

At non-zero T it is quite natural to expect non- 
vanishing thermal fluctuations of PC'^. Somewhat less 
trivial is the limit T — when the system approaches 
its (non-degenerate) ground state. In this limit no PC 
fluctuations occur only provided the current operator / 
commutes with the total Hamiltonian H of the system, 
otherwise fluctuations of PC do not vanish even provided 
the system remains exactly in its ground state^. For in- 
stance, PC does fluctuate down to T = provided the 
ring is coupled to some external dissipative environment. 
This situation is encountered, e.g., within the model^ 
where it was demonstrated that, while the average value 
of PC (/) in the ground state decreases with increasing 
coupling with the environment, PC fluctuations remain 
non-zero and even increase with the coupling constant. 

In the latter example interaction with the dissipative 
environment produces quantum decoherence which, in 
turn, causes suppression of PC (/) even at T = 0. On 
the other hand, ground state fluctuations of PC can oc- 
cur also in the absence of any decoherence. E.g., the 
operators I and H do not commute for a quantum par- 
ticle on a ring in the presence of an external periodic 
potential. This model describes, e.g., the properties of 
superconducting nanorings with quantum phase slips^. 
In this case PC fluctuations do not vanish even in the 
ground state* and, at the same time, quantum coher- 
ence remains fully preserved. As a result, the magnitude 
of PC fluctuations turns out to be sensitive to external 
magnetic flux*. This observation implies that quantum 
coherence and decoherence in meso- and nanorings can 
be probed by measuring the equilibrium current noise in 



such systems. The main goal of this paper is to the- 
oretically analyze such current noise in the presence of 
quantum decoherence produced by electron-electron in- 
teractions. 

Note that the existing theory of quantum decoher- 
ence of electrons in disordered conductors with electron- 
electron interactions?'- is rather complicated, merely be- 
cause of the Pauli principle which needs to be properly 
accounted for in such systems. At the same time, the ba- 
sic physics of this phenomenon can be explained already 
without unnecessary complications: It is due to the elec- 
tron interaction with the fluctuating quantum electro- 
magnetic fleld produced by other electrons moving in a 
disordered conductor. Hence, for the sake of physical 
transparency it is sometimes useful to employ a simpli- 
fled model which mimics all key features of the "real" 
problem of interacting electrons in a disordered conduc- 
tor except for the Pauli exclusion principle. This model^ 
deals with a quantum particle moving in a ring and inter- 
acting with some quantum dissipative environment. The 
latter could be, e.g., a bath of Caldeira-Leggett oscilla- 
tors or electrons in a disordered conductor. In the case 
of Caldeira-Leggett environment decoherence of a quan- 
tum particle on a ring was investigated with the aid of 
imaginarj*^ and real-time^^ approaches which yield sim- 
ilar results, i.e. exponential suppression of quantum co- 
herence down to r = at sufficiently large ring radii. 
In fact, the problem^ii^ is exactly equivalent to that of 
Coulomb blockade in a single electron box where expo- 
nential reduction of the effective charging energy at large 
conductances is also well established* ^^^iS,. The model of 
a particle in a diffusive electron gas was extensively used 
by different authors^ii^"— in order to investigate the ef- 
fect of interaction-induced decoherence on the average 
value of PC. Below we will employ the same model in 
order to study an interplay between PC fluctuations and 
quantum decoherence. 

The structure of our paper is as follows. In Sec. 2 we 
define our model and describe our basic formalism of the 
influence functional. In Sec. 3 we analyze PC fluctua- 
tions within the framework of the perturbation theory in 
the interaction. In Sec. 4 we go beyond the perturbation 
theory and evaluate equilibrium current-current correla- 
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tion functions in the limit of strong interactions. Our 
main conclusions are outlined in Sec. 5. Further techni- 
cal details of our calculation are specified in Appendices 
A and B. 



II. THE MODEL AND BASIC FORMALISM 

In our analysis we will employ a model^'''^'^ of a quan- 
tum particle with mass M on a Id ring of radius R pierced 
by magnetic flux ^x- This quantum particle interacts 
with a dissipative environment (bath) described by some 
collective variable V representing its degrees of freedom. 
The total Hamiltonian for this system reads 



H = 



+ Hcnv{V)+Hint{0,V), (1) 



where is the angle variable parameterizing the particle 
position on the ring, (f> = —i-^ is the angular momen- 
tum operator, = $/<i>o and $o is the flux quantum. 
The first term in Eq. ([T]) represents the particle kinetic 
energy, He^viV) is the Hamiltonian of the bath, and the 
term Hint{0, V) describes interaction between the particle 
and the bath. Note that within our model the interac- 
tion term involves only the coordinate (not momentum) 
operators. 

Let us define the current operator for a particle on a 
ring. It reads 



e ~ *e ~ - eld) — (b^ 

I = —e = —lH,e] = — — 

277 2tt^ ^ 2ttMR^ 



(2) 



Employing the Heisenberg representation of this operator 
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(3) 



one can define the current-current correlation function 
{I{t)I{0)). Evaluating the symmetrized version of this 
correlator in equilibrium one arrives at PC noise power— 
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did 



s(t) = -{i{t)m+mm) {if = j ^s^e-^-\ 

(4) 

which represents the central object in our subsequent 
analysis. 

Similarly to Refs. ilH we will model the environment 
by 3d diffusive electron gas with the inverse dielectric 
function 
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Dk'' 



e{u!, k) 



Ana 



(5) 



where a is the Drude conductivity of this gas, D = vpl/Z 
is the electron diffusion coefficient, w^? is Fermi velocity 
and I is the electron elastic mean free path. As usually, 
below we will assume disorder to be not very strong im- 
plying this mean free path to be much longer than the 



inverse electron Fermi momentum k'^ , i.e. kpl ^ 1. On 
the other hand, the electron mean free path should re- 
main much smaller than the ring radius I ^ R. We also 
note that Eq. ([S]) applies at frequencies lu lUc ^ vp/l. 

Fluctuating electrons cause fluctuations of the electric 
potential V in the system. Within Gaussian approxima- 
tion such fluctuations are described by the correlator 



(FFW = -coth|,Im^j^. 



(6) 



Interaction between the particle on a ring and fluctuating 
electrons in the environment is described by the standard 
Coulomb term 



Hint = eV, 



(7) 



where e denotes the particle charge. 

In order to evaluate the correlation function Q it is 
necessary to describe quantum dynamics of our system. 
For this purpose let us introduce the evolution operator 
U{t,t()) and define the density matrix operator 



(8) 



where pi is the initial density matrix. The kernel of the 
evolution operator can be expressed as a path integral 
over the angle variable. We have 

OO 

X ^ e2'^""'^-f/(6'i -h27rm,t;6li,0), (9) 



m— — OO 



where 



Ui9i,t;9[,0) = J VVVOe 



(10) 



e(o)=e; 



Sonv is the action of the environment and 




eV{vg{h),ti)\dti. (11) 



Here Ec — 1/{2MR^) and vg is the vector with compo- 
nents (i? cos 6',i? sin 6*). 

As we are interested in the dynamics of the particle 
rather than that of the bath, it is convenient to trace out 
fluctuating potential of the environment V . Making use 
of a standard simplifying assumption that at the initial 
time moment the total density matrix is factorized into 
the product of the equilibrium bath density matrix and 
some initial particle density matrix pj, we obtain 
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mi,m2 — — oo Q 
e-^(t)=ei+2irmi 0^(t)=02+27rm2 t 

X ve'^ei no^.e% (12) 

where p{6i,62;t) = {ei\p(t)\e2) and 

D / -»/(e"^'''(r8F(t')>t')-ey^(r„B(t'),t'))rft'\ 

-FK,0^]-(^e 0^ ^ey> >) \^ ^eM~^SR~Si), (13) 

is the influence functional^^ which depends on the angle variables on the forward and backward parts of the Keldysh 
contour, respectively 9^ and 6^ . Calculation of this influence functional amounts to averaging over the quantum 
variable V which is also defined on the Keldysh contour. This procedure^ can easily be adapted to our present situation 
of a particle on a ring where no Pauli exclusion principle needs to be taken into account, cf., e.gJ^. Introducing the 
new variables 9+ — {9^ + 9^)/2 and 9- = 9^ — 9^ , after the standard algebra (see Appendix A for further details) 
we obtain 

oo ' 

SR[9+,9^]^naY, 

a„n / dt'9+{t')sm{n9^{t')), (14) 

n=l { 

and 

oo * * 

Sj[9+,9_]^~27raY,a^J dt' J dt" ^.^^^.J^' ^„.. cos{n{9+{t') ~ g^(t-))) sin sin ^^g^, (15) 





where a — 3/{8kpP) is the effective coupling constant 
in our problem and a„ are the Fourier coefficients equal 
to a„ = (2/(7rr)) ln(r/n) for n < r = R/l > 1 and 
to zero a„ = otherwise. The weak disorder condition 
kpl 3> 1 obviously implies a ^ 1, i.e. the coupling 
constant always remains small within the applicability 
range of our model. 

We also note that the above influence functional re- 
duces to the Caldeira-Leggett one if we choose^ii^ a — 
r]B?/'K (where 77 is effective viscosity of the Caldeira- 
Leggett bath), ai — \ and a„ = for all n > 1. 



III. PERTURBATION THEORY 

Let us assume that both effective coupling constant 
a and ring radius R are sufficiently small and proceed 
within the perturbation theory in the interaction. Keep- 
ing only the first order correction a to the average 
value of PC, in the low temperature limit T <^ Eq one 
findsH 



for -1/2 < < 1/2. The resuh (HH) should be pe- 
riodically continued outside this interval of flux values. 
We observe that for a ^ 1 the first order perturbative 
correction to the average value of PC is negligibly small 
except in the immediate vicinity of half-integer flux val- 
ues (t)x = ±1/2, ±3/2, where the two lowest energy 
levels become very close to each other and the perturba- 
tion theory fails already in the first order. 

Turning now to PC noise let us recall that in the limit 
a — the current operator ^ commutes with the to- 
tal Hamiltonian. Hence, in the absence of interactions 
PC noise vanishes in the low temperature limit for the 
model in question^. On the other hand, in the presence 
of interactions PC noise remains non-zero and important 
contributions to the current-current correlator (j4]) are ex- 
pected to occur already at sufficiently small values of a 
and R. This limit will be studied perturbatively below 
in this section. 



A. Density matrix and Dyson equation 
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^Ig") In order to proceed it will be convenient for us to pass 

to the momentum representation. Performing the Fourier 
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transformation of the density matrix 



im9i~\-i7i92 



(17) 



and making use of the influence functional in the form 
(|A11|) (see Appendix A), we obtain 



p{mi,ni2;t) ^ ^ p.ijn^.m^) 

m'^ ,7712— — CXD 

X (k{mi,t; m[,0)K*{m2,t; m'^, 0) 



(18) 



where 



271- 



e(t)=ei 



e{o)=e[ 

Let us expand the exponent in Eq. (jl9p in series in 
and evaluate Gaussian path integrals in all terms of this 
expansion. In the zeroth order (all Vn — 0) we obtain 



K{m, t; m', 0) = Sm,m'Ko{m; t) 



(20) 



with Ko{z;t) = (,-iEc{z-cl>^ft ^ ^^^.^^ ^-^^^i vi ^ I and 
Vn^i yields the contribution 



Sm-n,m'KQ{m; t - ti)Ko{m - n; ti), 



(21) 



and so on. Representing each term in this expansion di- 
agrammatically, we can indicate the unperturbed propa- 
gator 5m,m'Ko{m; t) (I^Hl) by a solid line and observe that 
each insertion of e*"^*^*' (or e~*"^^*^) adds (or removes) 
the momentum n at a time t. Collecting all contribu- 
tions to all orders in a in Eq. (IT51) and averaging over 



Vn we arrive at the perturbation series for the density 
matrix which can be expressed in terms of Keldysh dia- 
grams consisting of two solid lines (implying forward and 
backward propagators) connected by dashed lines corre- 
sponding to the propagators IIJ^'"" {t — t') for the z/„-field. 
One of such diagrams is depicted in Fig. [T] and is similar 
to Keldysh diagrams encountered, e.g., in the Coulomb 
blockade problemSi. 

Let us denote the sum of all diagrams contributing 
to the evolution kernel for the particle density matrix 

as ^ml'l^!'{t), where lower (upper) indices correspond to 
forward (backward) lines. Provided upper and lower in- 
dices coincide, only lower indices will be indicated, i.e. 
^m'm' (^) = Um,m'{t)- The latter quantities can also be 
considered as elements of the matrix U{t). In what fol- 
lows we will also use the short-hand notation for the den- 
sity matrix V^^it) = p{mi,m2;t) and denote its diago- 
nal elements as Vm{t) or as a ket- vector \V{t)). In this 



m 



m-n 



m 



m-n 



m-n 



FIG. 1: Typical Keldysh diagram contributing to the density 
matrix. 



notations the evolution of the density matrix is expressed 
by means of the equation 



(22) 



Introducing the self-energy 5™^'™,^(i — t') in a standard 
manner as a sum of all irreducible diagrams we arrive at 
the following Dyson equation for the kernel of the evolu- 
tion operator 



U. 



xi^o(™i; t - h)K*im2;t ~ ^2)5™?;;:? (^2 - ii)^;™; (^i)' 
I 



(23) 



Suppose that our system was described by diagonal 
density matrix at some initial time which tends to —00. 
Then it should remain in the diagonal state at all later 
times as well. The evolution of such density matrix is 



determined by the equation 



d\m) 

dt 



dt'S{t-t')\V{t')). 



(24) 
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At long enough times \V{t)) tends to its equilibrium (and, 
hence, time-independent) value \T"^'^) which obeys the 
equation 



m-n 



m 



dt'S{t-t')\V'"i) = 0. 



(25) 



Performing the Fourier transformation of the self-energy 
Soj = j dte^'^*S{t) one can rewrite the above equa- 



tion as SqI'P'^'^) — 0. In equilibrium one obviously has 
-p^? = e--^'"/^/Z, where Z = ^e"-^"/'^ is the partition 

m 

function and _E,„ — Ecim — (jixY ■ 



B. Self-energy 

The evolution kernel for diagonal elements of the den- 
sity matrix can also be expressed via the self-energy as 



(26) 



Since all elements of the ket-vector \V{t)) are real, one 
can demonstrate that the self-energy remains purely real 
in the zero frequency limit. Introducing the matrices Ti^ 
and (which are non-singular at small frequencies) we 
can write 



(27) 



and, hence. 



a; + i(l + E„)-ir, 



-(1 + Sc..)"'- (28) 



Here we employ a simple approximation which amounts 
to neglecting S^^ and keeping only the leading in a cor- 
rection to r^. Then we obtain 



, -1^(0) ■ 
w -t- iTlj 



(29) 



where the real part of the self-energy is defined as 



u] \m+n.rn 



Traa\n\ I E„i+n — Em + IjJ 

^ Ve T - 1 



, (30) 



(0)1 



OO 



(0)1 



uj \ni,m / ^ \ uj Jm+n,m \^u} ^]m — n,m ) ■ (^l) 



The corresponding first order self-energy diagrams are 
depicted in Fig. [51 
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FIG. 2; First order self-energy diagrams 

C. Current-current correlator 

Let us identically rewrite the current noise power ^ 
in the form 

OO 

= 25R j dte'^'{I+{t)I+{Q)) - 27r5(c^)(/+)^ (32) 



where we defined /+ = (/^ + I^)/2. A typical diagram 
contributing to this expression is depicted in Fig. [31 Let 
us recall that in our problem the current and momen- 
tum operators coincide with each other up to a constant. 
Hence, inserting the current operator at a time t inside 
the diagram yields the factor eEcim — (j)x)/'i^^ where m 
is the particle momentum value at this time. In what fol- 
lows we will divide all diagrams into two different classes, 
those with equal momenta at both upper and lower at 
i = and all others. Summing up all diagrams in each of 
these two classes we arrive at two different contributions 
to the current-current correlator: 

(/+(t)/+(o)) = (£;|i(")w(t)i(°)|p^9) 

dt" [ dt'{E\I^^'^U{t~t')i^^^\t\t")\V^'>), (33) 



-OO 



where I*^"^ is a diagonal matrix with elements im]n — 
2Ec{m — 4>x)5m.n and I*^^^ (t, t') is defined as a sum of all 
irreducible diagrams containing i = 0. After the Fourier 
transformation we obtain 



^2TT5{uj){E\±^^^\V^'if. 



(34) 



This equation defines a formally exact expression for 
PC noise power. Proceeding perturbatively in a one can 
check that the contribution containing ii^' can be ne- 
glected as it turns out to be small as ^ a ^ 1 as com- 
pared to the terms with \ Then after some manipu- 
lations (see Appendix B for further details) we get 

= 2SR(^;iiw (u^- i(o)|7?-9). (35) 
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FIG. 3: Typical diagram contributing to the current-current 
correlator ([32}. 
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FIG. 4: PC noise power at T = ttq = 0.05, r = 5 and 
different flux values. Frequency uj and noise power Suj are 
normalized respectively by Ec and by e^Ec/{4:TT^). 



Note that at small lu the quantity tends to 

i\'P'^i)(^E\/{uj + iO), i.e. PC noise power is regular in 
the zero frequency limit. At non-zero uj the difference 
between these two terms is proportional to a, thus pro- 
viding nonvanishing PC noise at such frequencies in the 
presence of interactions. 
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FIG. 5: PC noise power at T = 0.05Ec and different flux 
values. Units and parameters are the same as in Fig. (4] 
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FIG. 6: PC noise power at (px ~ 0.15 and different tempera- 
tures. Units and parameters are the same as in Fig. |4l 



D. Results 



The perturbative expression for the current noise 
power psp was evaluated numerically at different tem- 
peratures and flux values. The results are presented in 
Figs. |3][S] One observes that the noise power strongly 
depends on the magnetic flux (p^- This property illus- 
trates the coherent nature of PC noised. PC noise grows 
with increasing and diverges as the flux approaches 
the point (jj^ = 0.5. This divergence has the same phys- 
ical origin as that in Eq. ((TB| . In this limit the dis- 
tance between the two lowest energy levels SE{(j)x) = 
Ec{l — 2 1 I) becomes small and the system undergoes 
rapid transitions between these energy states. As these 
levels correspond to different PC values such transitions, 
in turn, yield strong current fluctuations. 

It is important to observe that PC noise persists down 
to T — > 0. In this case S^j remains zero at frequencies 
smaller than the inter-level distance uj < 5E{(f>x) and 
becomes non-zero otherwise. We also note that in the 
lowest ('^ a) order of the perturbation theory zero tem- 



perature PC noise vanishes at (jj^ — 0. In general, how- 
ever, this feature does not hold, as it can be observed, 
e.g., from the expression for PC noise power in terms of 
the exact eigenstates of the total Hamiltonian^. Non-zero 
PC noise a.t 4>x = will also be demonstrated in the next 
section where we employ non-perturbative quasiclassical 
analysis of the problem. 

At non-zero T there appears additional zero frequency 
noise power peak. This peak grows rapidly with in- 
creasing temperature and eventually assimilates all other 
peaks. As a result, at sufficiently high temperatures 
only a wide hump remains, and PC noise becomes flux- 
independent, i.e incoherent. The dependence of zero fre- 
quency PC noise on temperature is illustrated in Fig. [T] 
Interestingly, at low enough T this dependence turns out 
to be non-monotonous, while in the high temperature 
limit it approaches the linear dependence Sq oc T, as it 
will also be demonstrated in the next section. 
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FIG. 7: Zero frequency PC noise power So as a function of 
temperature for (l>x = 0.15, na — 0.05 and r = 10. Tempera- 
ture and PC noise power are normalized respectively by Ec 
and by e'^Ec/i^TT^). 



IV. NON-PERTURBATIVE ANALYSIS 

Let us now turn to the limit of strong interactions 
in which case the effect of a dissipative environment on 
the particle motion becomes large substantially reducing 
fluctuations of the angle variable 0. It is important to 
stress that for the model under consideration this situa- 
tion can be realized even at small values of the effective 
coupling constant a <C 1 provided the ring radius be- 



comes sufficiently largei^, i.e. 



(36) 



In this limit and provided temperature is not too low it 
suffices to employ the semiclassical approximation and to 
expand the effective action (|14l) . ([T5|) up to quadratic in 
6'_ terms. As usuallj«22ri2ij the resulting effective action 
can be exactly rewritten in terms of the quasiclassical 
Langevin equation for the " center-of-mass" variable 9+. 
For the model studied here this Langevin equation takes 
the form 

—e+{t) - l9+{t) = ^(e„WcosM+(i)) 



2Ec 



-Xn{t)sm{ne+{t))), (37) 



where we defined 

oo 

7 = 27ra a„n^ = Anar^ 



(38) 



and introduced Gaussian stochastic fields £,n{t) with the 
correlators 



sinh^(7rT(t-t'))' 



(39) 



(40) 



In the high temperature limit these correlators reduce to 
those describing the white noise 

(e.(i)em(t'))«,A = 2Srn,„7:aann^TSit - t') (41) 

and the Langevin equation can be solved exactly. As a 
result, we arrive at the high temperature noise power 



(42) 



At lower temperatures the white noise approximation 
(|1T|) becomes inaccurate and Eqs. ([Mt . PO)) should be 
employed. In this case noise terms in the Langevin equa- 
tion can be treated perturbatively24. Keeping only the 
zeroth and the first order contributions one gets the so- 
lution of Eq. (|57)) in the form 



(43) 



where 9^^^ is an arbitrary (and physically irrelevant) con- 
stant and o'"^^ [t] obeys the equation 



2Ec 



5(1) 



n=l 



(44) 



Resolving this equation one immediately arrives at the 
noise power in the form 



c 



27r2(a;2 



,„-wcoth — , (45) 
i-fEcf) 2r' ' 



which again reduces to Eq. (I42p in the high temperature 
limit T ^ uj. Note that for w ^ jEc the parameter Ec 
drops out and the noise power becomes 



TT^coth— , 



(46) 



i.e. in this case S^j on 1/a. For a; — > this expression 
further reduces to So oc T/j. The noise power Scu ([45|l is 
also depicted in Fig. |S]at different values of T. 

Comparing Eqs. (|42|) . (|45|) with perturbative in the 
interaction results for the noise power derived in the pre- 
vious section we observe a striking difference between 
them: While in the weak interaction limit PC noise is 
sensitive to the externally applied magnetic flux (j)^., in 
the opposite limit of strong interactions the noise power 
S'lj turns out to be essentially independent on The 
latter observation implies that in the non-perturbative 
limit (j36p quantum coherence of the particle is suppressed 
by strong interactions with the dissipative environment. 
This conclusion is fully consistent with earlier results^^ 
derived for PC (/) in the limit (|36|) . 

Note that Eq. P5)) defines only the dominating con- 
tribution to the noise power. In addition there also exist 
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FIG. 8: Noise power at different temperatures for ttq = 0.05 
and r — 25. Units are the same as in Fig. 4. 



to do with the fact that the angle variable 9 is compact 
(i.e. defined on a ring). On the other hand, the Langevin 
equation method employed above effectively " decompact- 
ifies" our problem, thus being able to capture only the 
(/)2:-independent contributions to S^j- In order to esti- 
mate the leading 0a;-dependent correction to Eq. (|45|) we 
will make use of the approach initially developed for the 
problem of weak Coulomb blockade in metallic quantum 
dots2^"— . This approach establishes the relation between 
the density matrices and expectation values evaluated for 
the problems described by the same Hamiltonian but re- 
spectively compact and non-compact variables. With the 
aid o&rU for the expectation value of the current oper- 
ator one finds 



{!) = 



(47) 



small corrections to this result which do show the depen- 
dence on the external flux (px piercing the ring. Tech- 
nically, the existence of this flux-dependent terms has 



Here p„p stands for the reduced equilibrium density ma- 
trix for a particle described by non-compact (i.e. de- 
fined on a straight line) variable 9. Analogously the noise 
power is given by the autocorrelation function 



J 



S{t) 



{eEcY En tr((#„p + Pnp^Mpit, 0)<^e2-^(^-^^)[/„p(t, 0)) 



27r2 



^^tr(e2-w(^-^^)p„p) 
I 



(48) 



where the evolution operator C/„p(i,0) is again defined 
for a non-compact variable 9. With the aid of the path 
integrals one can rewrite the above equations respectively 

as 



In the semiclassical limit averaging in these equations is 
conveniently performed within the above Langevin equa- 
tion technique. With the aid of Eqs. (|43| and (|44|) one 
easily finds 



and 

Sit) 



277i7V(e+(t)/(2Ec)-0x) 



2tt 



^^(^+(O^+(O)e2-^^(^+(*)/(2i^c)-0.)) 



47r2 



^j^(e2"W(e+(*)/(2Bc)-0x)) 



(49) 



(50) 



_ ieK{Q) ^^jVe-2^»JV0.-^^jv^g(o)/(2g^) 



2Ec 



E 



N' 



-2iriNcl>^—K^N^K{0)/{2E^) 



and 



4.2 



K{t)( TT^KjO) Eat jV2e-2"^0.-^^^^^(o)/(2i?^) \ 



(52) 



where we introduced the correlator K{t) — We also obtain 
{9''^\t)9^_i^\{))) which Fourier transform equals to 



oo 

K^= J dte"^*Kit) = — 



27^c ,^eoth^. (53) 



^W = 47i^g./2^ ,2^(>^)2 ^oth^. (54) 



2T' 



2T 



Logarithmic divergence contained in this integral can eas- 
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ily be cured if we recall that (a) our diffusive electron 
gas model ([5]) is applicable only at frequencies w ^ Wc ^ 
vp/l (hence, the integral in Eq. (|54|) should be cut at 
uj ~ Wc) and (b) our Langevin equation approach be- 
comes insufficient in the low temperature limit where it 
should be supplemented by other techniques. As a result 
of these considerations we may write 

K{Q) = C + 2EcT 

where ip{z) is the digamma function and the constant C 
effectively accounts for the low temperature behavior of 
our system. The value of this constant can be determined 
if we compare the expression for PC ([ST|) derived here 
with the results of the instanton analysis^^. Since in the 
limit (p6| we have K(0) ^ E^, it suffices to keep only 
the terms with TV = 0, ±1 in Eqs. and Then 
comparing Eq. (jSip with the resulti^ (/) oc exp^—Anar) 
one may identify C as 

C^^. (56) 

TT 

Finally, for the current noise power we obtain 
_ e'^-fEl,ucoth{uj/2T) 

X (^1 ^f^e'"^ cos(2^</..) j . (57) 

As it was already anticipated, in the limit of strong in- 
teractions ([55]) the coherent (flux-dependent) contribu- 
tion in Eq. (j57p just represents a small correction to the 
main incoherent term ()45|) . In this respect an accurate 
evaluation of this small correction may even be consid- 
ered as exceeding, it suffices to demonstrate that this 
(/f^:- dependent correction remains small in the limit p6p. 
We also note that exponential dependence of PC on the 
ring radius (/) oc exp(— 47rar) applies down to temper- 
atures T ^ Ec / {4:Trar), whereas at even lower T ^ 
it crosses over to a weaker (power law) dependence (cf. 
Figs. 1 and 2 in Ref. [l3) but remains strongly suppressed 
in the limit ar ^ 1. Accordingly, one can expect that 
in the same limit the flux-dependent correction to the 
incoherent noise term (j45p remains small down to T = 0, 
though at r ^ Ec / {4:'!Tar) it may deviate from the form 
(j57p . Unfortunately, quantitative non-perturbative anal- 
ysis of the exact zero temperature limit appears difficult 
since neither Langevin equation approach nor instanton 
analysiai^ can be trusted in this limit. 

V. CONCLUSIONS 

In this paper we analyzed fluctuations of persistent 
current produced by a charged quantum particle moving 



in a ring and interacting with an environment formed by 
3d diffusive electron gas. Specifically, we restricted our 
attention to PC noise and evaluated symmetric current- 
current correlation function in two different limits of 
weak and strong interactions. Note that although within 
our model the effective coupling constant a describing 
Coulomb interaction between the particle and the bath 
always remains small, a <C 1, interactions can be treated 
perturbatively only for sufficiently small values of the ring 
radius R, while for larger R (j36p non-perturbative anal- 
ysis of interaction effects becomes unavoidable. 

In the absence of interactions within our model PC 
fluctuates only at non-zero T and no such fluctuations 
could occur provided the system remains in its ground 
state at T = 0. In the presence of interactions the cur- 
rent operator does not anymore commute with the total 
Hamiltonian of the system and fluctuations of PC gen- 
erally persist down to zero temperature*. In the pertur- 
bative regime of weak interactions and at sufficiently low 
T quantum coherence of the particle remains preserved, 
PC noise is coherent and, hence, the noise power can 
be tuned by external magnetic flux (jj^- In contrast, in 
the limit of strong interactions (j36p fluctuations in the 
electronic bath strongly suppress quantum coherence of 
the particle down to T = 0. In this case the average 
value of PC (/) gets strongly suppressed as welU^, while 
the current noise, on the contrary, does not vanish and 
becomes practically flux- independent. In other words, in 
this regime fluctuations in the environment induce inco- 
herent background current noise in the ring which per- 
sists even at zero flux cj)x — when the average PC is 
absent (/) = 0. 

We also point out that, while in the perturbative limit 
PC noise power Suj tends to increase with the coupling 
constant a, in the non-perturbative regime the depen- 
dence of Suj on a becomes more complicated, cf. Eqs. 
(|i5)l and ([5^ . In particular, at sufficiently low frequen- 
cies we find Suj oc l/{aR'^), i.e. in this regime the noise 
power decreases with increasing both a and the ring ra- 
dius R. On the other hand, the average PC value de- 
creases even much stronger and, hence, the ratio S^/{I) 
increases with increasing a and R. 

Perhaps the most important result of this paper is the 
prediction of (i) coherent flux-dependent fluctuations of 
persistent current in sufficiently small rings and (ii) in- 
coherent flux-independent current noise in larger rings. 
Thus, quantum coherence and its suppression by inter- 
actions in meso- and nanorings can be experimentally 
investigated not only by detecting the average PC value 
(which can happen to be very small) but also by measur- 
ing PC noise and its dependence on the external magnetic 
flux. We believe it would be interesting to perform such 
experiments in the near future. 
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Appendix A: Influence Functional 

Let us present the derivation of the influence functional defined by Eqs. (fT3)) . (IT4l) and (|15p . In order to perform 
Gaussian averaging over fluctuating electric potential it is sufficient to define only the second order voltage correlators. 
Introducing the variables V'^ = {V^ + V-^)/2 and V" = — we can express these correlators in terms of the 
Green functions 

{V+{re{t),t)V+{r^{t'),t')) = zG^(2i?| sin((^^ - x)/2)|, i - t'), (Al) 
{V-{rg{t),t)V+{r^{t'),t')) = zG^{2R\sm{{e-x)/2)\,t~t'), (A2) 

{V+irgit),t)V-ir^it'),t')) = zG^(2i?|sin(((?-x)/2)|,t-i'), (A3) 
{V-{ro{t),t)V-{r^{t'),t')) = 0, (A4) 

where the last equation is a direct consequence of causality. Here G^, and G^ are respectively retarded, advanced 
and Keldysh Green functions related to the dielectric function e(k, of the environment as follows 

G^(k,c.) = — -— G^(k,w) = (G«(k,c.))*. (A5) 

k^e(k, w) 

2G^'(k,L.) = coth^ (G^(k,w) -G^(k,c.)) . (A6) 
Combining the above expressions with Eq. ([S|), in the case of a diffusive metal we obtain 

G\X,t- 1') = , 5'{t - 1') G^iX,t- 1') = , , . (A7) 

Gaussian averaging over the fields can now easily be performed, cf., e.g.^. As a result we arrive at Eq. (|T3)) . The 
expression for the imaginary part of the action in this equation reads 

t t 



+C{e^{t') - o^{t")) - c{o^{t') - e^{t")) ~ cio^it') ~ 0^{t"))) , (A8) 

where (^{x) = [4(i?/Z)^ sin^(x/2) + 1]^^/^. Expanding in the Fourier series 

C{x) = ao - ^ a„ sin^ (^^"j = ^ + - ^ a„ cos(na;) (A9) 

71—1 n—1 

with an — (2/(7rr)) ln(r/n) for n < r and an — ^ for n > r and using the identity 

6 6 

cos(6'f - e^) + cos(6lf - ) - cos(6if - 0^) - cos(6if - 6^) = 4cos(6'i+ - 02+) sin -j- sin ^ (AlO) 

one arrives at Eq. (fT5|) . Eq. (fT4|) is recovered in a similar manner. Let us also note that the Caldeira-Leggett 
environment is described by the function (^{x) = — 2sin^(.x/2) = cos(a:;) — 1 which should be employed in that case 
instead of Eq. (K9\ . 

Let us also rewrite our infiuence functional in a somewhat different form, which can be conveniently used in our 
perturbative calculations. Employing the definition (—1)''^ = 1 and ( — 1)^ = —1 we obtain 

J-[6i^,0^] = / e "=i« ) (All) 

Here Vn is Gaussian stochastic complex variable described by the correlator 

(-r(^)-™'(t')) = ( Dlit t') + t^D^^it' -t) + tl^Df^it - t')] . (A12) 
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where 



Then we obtain 



Dil it - t') = TTaan5'{t - t'), D'^ {t - t') 



-iiraan 



ttT' 



(A13) 



E E f dt j dt'ni-' {t-t')cos(n{0''{t)-e' (t'))) 

^=1 (T,cr'^F,B 



With U-'-'it - t') = -[-lY+-\yl*{t)v-^{t')). 



(A14) 



r 



Appendix B: Operations with singular matrices 

^ The evolution kernel for the diagonal density matrix 
can be expressed via the self-energy by means of the 
following equation 



(Bl) 



From the identity {E\Su = we conclude that the ma- 
trix has zero eigenvalue with the left eigenvector {E\. 
Hence, there also exists the right eigenvector \uj) with the 
same (i.e. zero) eigenvalue, Su:\oj} = 0. Employing the 
normalization condition = 1 and introducing the 

projector jC^^ = \uj){E\ we can verify the identity 



i£,jCuj \ i 



(B2) 



Indeed, expanding the above expression at small frequen- 
cies we obtain 



Observing that at zero frequency the vector \lu) just co- 
incides with the equilibrium distribution ^function, i.e. 
|0) = jP'^'^) and making use of equations <So = —To and 
<Sq = iSo we get 



+ - 



i\P^''){E\ 

uj + iO 
1 



^\v^i){E\ + ro 



(l-(l + Eo)r'^)(£;|) (B4) 



which holds^ for any value ^. This identity irnplies that 

the matrix is singular at small frequencies, oc -^j^- for w — > and any value of ^. 
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